Non-Gaussian features of primordial magnetic fields in power-law inflation 
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We show that a conformal-invariance violating coupUng of the inflaton to electromagnetism pro- 
duces a cross correlation between curvature fluctuations and a spectrum of primordial magnetic 
fields. According to this model, in the case of power-law inflation, a primordial magnetic field is 
generated with a nearly fiat power spectrum and rms amplitude ranging from nG to pG. We study 
the cross correlation, a three-point function of the curvature perturbation and two powers of the 
magnetic field, in real and momentum space. The cross-correlation coefficient, a dimensionless ratio 
of the three-point function with the curvature perturbation and magnetic field power spectra, can 
be several orders of magnitude larger than expected as based on the amplitude of scalar metric 
perturbations from inflation. In momentum space, the cross-correlation peaks for flattened triangle 
configurations, and is three orders of magnitude larger than the squeezed triangle configuration. 
These results suggest likely methods for distinguishing the observational signatures of the model. 
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I. INTRODUCTION 

The simplest models of inflation have proven remarkably successful at solving an array of cosmological problems 
whilst bringing to order a vast catalog of observations. Hovi^ever, little is known about the Lagrangian during infla- 
tion. Non-gaussian correlations of cosmological observables carry information about the details of the physics during 
inflation. In this article, we study what sort of dominant non-gaussian features should be expected in the statistics 
of primordial magnetic fields arising from a coupling between the inflaton and electrodynamics. 

The Lagrangian we examine is motivated by the problem of galactic magnetism. It is observed that galaxies 
are permeated with magnetic fields of order 100 to 1 fj,G. This can be seen in great detail in the recent study of 
magnetic fields in M51 by Ref. [l[, for example. The current understanding of the formation and evolution of these 
fields is that a dynamo mechanism, combining galactic rotation with helical turbulence, amplifies a seed magnetic 
field with strength larger than 10"^" G at the time of the galactic disk formation 0-3- The seed field must have 
been more or less homogeneous on the scale of the galactic disk when galaxies started to form, implying a comoving 
coherence scale of around a Mpc. Such a large scale has been inside the horizon since a redshift z ~ 10®, however 
the physics from the time of Big Bang Nucleosynthesis at z ~ 10^ to the present is well-established and, in the 
standard cosmology, no large-scale magnetic field of significant strength is known to have been generated inside the 
horizon since then. In fact, the largest seed field produced from the nonlinear evolution of subhorizon primordial 
density perturbations, while remaining consistent with the physics of recombination, is only 10^^^ G on Mpc scales 
at the present time iSj, which is far too small to seed the dynamo mechanism. Logically, then, we may ask whether 
a superhorizon magnetic field produced during inflation could have provided the seed field. However, quantum 
electrodynamics is conformally invariant, so that the amplification mechanism that produces density pertubations 
during inflation leaves an uninterestingly small magnetic field energy density. To explain such a seed field, new 
physics would be required during infiation to break conformal invariance [6|. 

In a previous article [3], a toy model for such new physics was investigated in which a spectator field in de Sitter 
spacetime couples to electromagnetism. In the absence of gravity, it was shown that a scale-free magnetic field 
spectrum with rms amplitude ^ 0.1 nG at Mpc scales may be achieved. The central result of that investigation was 
the cross correlation between the scalar field and magnetic fields, showing that the dimensionless amplitude, measured 
in units of the power spectra, can grow as large as ~ 500Hi/M, where Hi is the inflationary Hubble parameter and 
M is the effective mass scale of the coupling. 

In this article, we consider a more realistic scenario, coupling the inflaton of power-law inflation to electromagnetism 
with a dilaton-type interaction. Many of the results from the earlier investigation will be shown to carry over, including 
the result of a nearly scale-invariant magnetic fleld spectrum with rms amplitude ranging from 10"'^ — 1 nG at Mpc 
scales. However the present work concentrates on the issue of the non-gaussian corrections to the statistics of the cosmic 
magnetic field [8[ . These are contained in terms such as the one-loop correction to the magnetic field power spectrum. 
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four-point function, or the cross correlation with metric perturbations. The first two involve at least two insertions 
of the interaction Hamiltonian in the in-in formalism perturbative series, hence are expected to be subdominant in 
comparison to the latter. Since tensor perturbations are harder to detect than scalar metric perturbations, the most 
important non-gaussianity in magnetic field spectra arises in the cross correlation between the curvature fluctuation 
7?, and the magnetic field energy density, in the form (JZB"^) which we evaluate in this paper. This cross correlation, 
which we naively expect to have amplitude 10^^ to match the amplitude of scalar metric perturbations, may be probed 
directly in experiments that are sensitive to the statistics of the large-scale, initial conditions in the mass density and 
magnetic fields, such as cosmic microwave background temperature-polarization correlations and the correlation of 
all-sky Faraday rotation measures with galactic density 9| . 

The outline of the paper is as follows. In Sec. |IT]we introduce our model. In Sec. IIIII we present the calculation of 
the cross correlation, and we study its properties in Sec. IIVI Wc conclude in Sec. |Vl In regard to notation, a prime 
indicates a derivative with respect to conformal time, ' = d/drj, an overdot indicates a derivative with respect to 
cosmic time, ' = d/dt, and dt/dr] — a where a is the expansion scale factor. When necessary, an asterisk subscript 
denotes a quantity evaluated outside the horizon. 

II. RELIC MAGNETIC FIELDS FROM INFLATON 

A. Amplification mechanisms 

In an effective field theory of inflation, it is possible to couple the inflaton ip to electrodynamics at the lowest 
dimensionality in general form with an action 



S= d 
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M^ 1 1 

Idtt 2 4 



(1) 



where VF is a dimensionless scalar built out of Lp and its derivatives. We envisage an ultraviolet completion of the 
effective Lagrangian in Eq. ([1]) where integrating out heavy fields leads to the infiaton potential V{tf) and the coupling 
W{ip). In the very low energy limit, we require W to be constant, so that quantum fiuctuations of the electromagnetic 
field are conformally invariant, and there is no coupling between electrodynamics and the inflaton. However, conformal 
invariance is broken when W evolves in time or varies in space. 

In a previous work, we showed that by driving W with a spectator field during infiation, then the evolution W (x a^ 
produces a relic cosmological magnetic field as large as 10^^° G on Mpc scales at the present time |7|. A similar 
approach was previously studied in Refs. [10l - ll2l |. Our function M^ is a special case of the function denoted P in the 
general case studies of Refs. |13l.ll4|. Our plan in this work is to show that once an amplification mechanism has been 
selected and the power spectrum of magnetic fields is used to fix the parameters of the Lagrangian in Eq. ([1} , one 
can go further and predict the structure of the cross correlation of scalar metric fluctuations and magnetic fields on 
cosmological scales. 

As a concrete effective theory, we consider slow-roll inflation with the potential [15l - [l8| 



V{Lp)^M^eyiY>{-Lp/k) (2) 

where M and A are two mass scales. In a theory where exp(iy9) appears in the Lagrangian, it is natural to take as the 
coupling function 

I^(^)=exp(^) (3) 

where e = —H/H^ ~ l/(167rGA^) = Mpj/(167rA^), n is a constant parameter, and we indicate the Planck mass as 
-Wpi — ^l\G. The appeal of the single-field inflaton potential of Eq. ([2]) is that it admits a closed form solution for 
the scale factor and the evolution of the fields in time without any further approximation in the size of e, namely 
H = 1/et, a{t) = ai{t/tiY''^^ or 0(77) = ai{rii/riY'^^~'^^ in terms of conformal time, where the subscript / indicates 
values at the end of infiation, and 



^W=Al°g(7:^-wF)- (4) 



Therefore, the coupling function is 

W{^) = Wi{alaif- (5) 



where Wi is the couphng at the end of inflation. 

For arbitrary inflation models, it is always possible to start with a generally covariant gravitational form like Eq. ([3]) 
and arrive at the specific slow-roll behavior of Eq. ([5]). This is because the scale factor for generic single field slow-roll 
evolves according to 



a{t) ~ exp 
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(6) 



and the fiatness condition requires that V{(p)/V'{(p) is approximately a constant as a function of ip, hence a{t) will 
be a power of exp{(p/v). Therefore, W will be proportional to a if it is a power of exp{(p/v). Then, as reheating takes 
place, the inflaton fleld value and its derivatives approach zero, so that VF — > 1. Consequently, there is no modiflcation 
of electrodynamics at later stages (or lower energies). 



B. Quantum fluctuations of the electromagnetic field 

The action Eq. (HI) of the free electromagnetic fleld is 



'4^/ — -jxr/.^\^ I? J7P.V - / J», J3^ T/t/, ^^^^^ /„/'^,^^2'^ ( 1^2 _ 1^2 



d^x^^W{v)-F^,.F^''' - / dj^d'x Wi {a{jj)/a{j^i)) 
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(7) 



where above and hereafter, Latin indices z, j, ... indicate components of the comoving Cartesian coordinate system. 
The quantum field theory will avoid ambiguities if we choose a complete gauge specification, thus we proceed in the 
Coulomb gauge defined by diAi = in these coordinates. Gauss' law is then satisfied with Aq = 0, and the action 
further simplifies to 



Sa ^ jdijd'x Wi {a{ij)/a{7jj)f- Q^f - ^{d^A.^^ 



(8) 



The canonical field conjugate to Ai is ni(x, rj) = 6S/SA'^{'k, rf) — Wi{a/ai)'^'^A'^{-x., rf), thus the canonical commutation 
relation [Ai{'x.,rj),Iij{y,r])\ = i5ij5^{x. — y) is here equivalent to 



[A,(x,7y),A;.(y,,7)] = ^{ai/af-5.,5H^-y) 



We perform the usual expansion 

A,(x,77)=^ 
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where /3(k, u) and /3(k, cr)^ are the annihilation and creation operators, which satisfy [/3(k, a),/3(k',(T')^] = 
(27r)'^(5CTo-'i5'^(k — k'). Moreover, the polarization vectors e,^"^ (k) satisfy Eq. ([9]) if 

^ei^)(k)e5.^)*(k) = %-a,. (11) 



Finally, 



VkivWkiv) ~ vlivWkiv) = -^("-r/")^" 



gives the Wronskian condition on the mode function Vk {rj) ■ 

The equation of motion for the gauge field mode function following from Eq. 

v'l + 2n—vi + k'^Vk = . 



is given by 



(12) 



(13) 



For successful inflation, all relevant cosmological modes must have been deep inside the horizon at an early enough 
time: 
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(14) 



which is consistent with the Wronskian condition, Eq. p2)) . and fixes the sohition of Eq. p3|) we seek: 

'"'^'^^ = (l^) ^'"""^'^"^Hv/ViTHi^H-kv) (15) 

= Vk*Ua{-kr]) (16) 

-"(-) - ^^^'^(-) (17) 

where a — 1/2 + n/{l — e). The constant coefficient Vk* is the value of the mode function outside the horizon, 
Vkifj) = Wfc for —krj <^ 1 where 

_ . 2"-ir(a)(^fc,,)V^-" ,,^,,,, 

Thus, Mpc modes of the electromagnetic field are stretched beyond the Hubble scale during inflation, where they are 
frozen until re-entering in the present epoch. The result of Eq. (J18p will prove useful in determining the size of the 
magnetic field at the present time. 

C. Magnetic field power spectrum 

The magnetic field two-point correction function is given by 



1 /. d^ d 



2 



(5.(x,,)S'(y,,)) = -,[S.,^^;jg^~g-g-]{Aix,V)My^V^ 



(27r)3 
PBik) = 2%^ = l ''^tl^;/\ -k,fH^^\-k,)H^^\^k,). (19) 

The above spectrum is identical to the result obtained in our previous work, Ref. [7|, upon taking the limit e — > 0. 
The power spectrum for Fourier modes outside the horizon at the end of inflation, when \krji\ ^ 1, is given by 

A|(fc)^^P5 = U \^^\J {l-efHf{-krj,r^ (20) 



27r2 \ttJ 2^bWi 

where ub = 4 — 2n/(l — e) for n > 0. In the second line above, we express the power spectrum at the end of inflation 
in terms of the scale factor and Hubble parameter at horizon cross, namely k/a^ — H^. We can see that ns = 0, 
or n = 2(1 — e), gives a flat spectrum. (A flat spectrum can also be achieved for n = —3(1 — e) but this solution is 
not viable, as discussed in the next section.) After inflation, when W{ip) becomes a constant, conformal invariance is 
restored and the magnetic field amplitude will decay as a~^. Under the assumption of instantaneous reheating, the 
spectrum today will be {aj/ao)^ times the above result. 

To obtain a numerical value of the power spectrum, we start by fixing e ~ 0.016 and H^ ~ 10^^ GeV for Mpc modes, 
both values of which are consistent with WMAP. Under the assumption that the Hubble parameter does not change 
from its horizon cross value by many orders of magnitude during reheating, one may write 3H^ /SttG — n^giTj/SO 
where gj and T/ are the effective degrees of freedom and temperature of the radiation fluid at the end of inflation, 
and we assume adiabatic evolution thereafter. For this simplistic model, we use gi '^ 10^ to take into account possible 
degrees of freedom beyond the Standard Model. As a consequence, T/ ~ 3 x lO^^GeV and iJ/ ~ 5 x 10^^ GeV, 
corresponding to a redshift for the end of inflation of zj ~ 10^^. (Note that in Ref. \T\ we used Hj ~ 10^"* GeV 
and zi = 10^^.) With the aid of the conversion factor (Gauss)2/87r = 1.91 x IG^'^'GeV , the energy density per 
logarithmic frequency interval of magnetic fields at the present time is 



dlnA;' ' ^' ' \7: J 2"bW^/ ' ' * \ai J \ao J \a^H. 

{{b-nB)l2f ( k 
M^/r(5/2)2 l^Mpc 



10^22.8-22.5». r((5 ^i3)/2)V _fc I ^2 (23) 



If Ub = and Wj ~ 1 then the field strength is roughly 10~^^ G on Mpc scales, whilst a slightly blue tilt as great 
as Ub ~ —0.2 yields -Brms — 10~^ G, which in both cases may be sufficient to explain the observed astrophysical and 
cosmological magnetic fields. 



D. Backreaction 



Self-consistency of this scheme for the generation of a primordial magnetic field requires that the amplification 
mechanism does not generate an energy density of gauge fields as large as the inflaton energy density and thereby 
spoil inflation. This is a real concern, since the n = —3 case, while corresponding to ub = 0, leads to an overproduction 
of electric field energy density. This issue was first analyzed in detail by [ij] . As shown in our earlier work [7] , the 
energy density of the electromagnetic field as measured in the cosmic rest frame is 



W{r]) 
W{r]) 



PEB = ,^,{A[^, + {d,A,mAj)) (24) 



2.MV)'J '^''^'^d^^l' + '^'l^^l') (25) 

nEB - \GHJ{l~ef j'^_^^dxx' {\H[^\x)\^ + \h'^^\{x)\^) (26) 

- rFf2 f(^(l), for |n|<2 

where Ni is the number of e- foldings of inflation. Using Ni = 70 as a fiducial value, we derive the limit \n\ < 2.1. 
Using nB — ^— 2n/{l — e) for n > 0, this translates into the bound —0.2 < ns < 4 within which ^Ieb < 0.01, such 
that the energy density of quantum electrodynamics amplified by the inflaton will be small. This leaves the ub — 
case as a potentially viable scenario for the creation of a scale-free spectrum of primordial magnetic flelds. 

E. Charge energy density contribution 

There is a further issue to be addressed in a successful model of primordial magnetogenesis, which was raised in 
Ref. [IJI regarding the effect of the electromagnetic current. To illustrate the issue, suppose that we consider the 
charged fields at the relevant scale H^, to be complex scalars ct„: 



S = 



-\w{ip)F^^''F^, - J^A^ -Y^g^^'-d^ald^an 



(28) 



and J^ is the gauge-invariant covariantly-conserved current associated with the matter fields cr„. Based on the 
discussion of the previous section, Sec. IIIDl we would like to take n ^ 2 which means W^ ex a*. If the ratio of 
Coulomb energy to electromagnetic kinetic energy is of order one at the end of inflation then it is of order e^^' at 
early times, where Nj is the number of e-foldings, therefore breaking down perturbation theory. In other words, 
the electromagnetic coupling is proportionately large and therefore strongly coupled at the beginning of inflation. 
This would appear to be fatal for this model, since we could no longer justify a free-fleld, perturbative treatment of 
electromagnetism. Equivalently, if we start with W of order one at the beginning of inflation in Eq. (|28p . then at the 
end of inflation the kinetic energy of the gauge fields would acquire a factor e^-^' . If at this point we redefine the 
gauge field to canonical normalization, by scaling v4^ by a factor of e~^^' , then the classical electric charge becomes 
e~^'^'eo, where cq is the electric charge at the beginning of inflation. Since one expects iV/ ~ 65 — 70 at least, the 
electric charge at the end of inflation is at least 130 orders of magnitude smaller than one. This is to be taken as the 
charge at the scale of reheating. If the beta function of the electric charge at these scales is negative, then its value 
may be brought to larger values at zero energy, but this would require an amplification of a similar number of orders 
of magnitude. 

The problem of the strong coupling can be fixed if the effective field theory at the relevant scales breaks gauge 
invariance as follows: 



S = 



-\w{ip)F^^'^F,,, ~ W^ipyj'^A^ -Y^g'^'^d^^d^o^ 



(29) 



With the convention W = Wi{a/aj)'^'^ and Wj ^ 1, the action at the end of inflation (a — aj) reads more expHcitly 



S^ I d'x dr, -lr,''Prj''^F^,Fp, " I] «? (ieo<9^(7„ - H.c.) 77^" A, - ^ ajrj^''d^,a*d,a 



(30) 



where eg is the electric charge during inflation. In this case, ignoring any running of the electric charge from reheating 
to zero momentum transfer, one can simply take cq to be the renormalized charge at zero momentum. There is no 
longer any problem at early times since both the current contribution and the kinetic energy of gauge fields go to zero 
as a —> at the same rate. As pointed out in Ref. [19\ this would not succeed if gauge invariance is preserved. 

In this work we will concentrate on the possible non-gaussian signature of the amplification mechanism and so we do 
not provide any rationale for writing a Lagrangian like Eq. (l29l) , which would have to arise from the UV completion. 
However, we point out that in theories with extra dimensions, gauge invariance can be violated in the 3+1 brane in a 
controllable way at high energies as a consequence of charges leaking to the extra dimensions, even if gauge invariance 
is maintained in the full space-time. Consider e.g. the discussion of this issue in Ref. [20|]. In this latter reference, it 
is shown how charged particles leaking to the extra dimensions produce electromagnetic waves in the 3+1 brane. It 
would be interesting to analyze whether this effect in brane-world models may provide a physical mechanism for the 
amplification factor W. 

F. Curvature perturbations 

Quantum fluctuations of the inflaton field create inhomogeneities that may be treated in terms of perturbations 
of the spatial curvature, TZ. A suitable set of the coordinates of the perturbed Robertson- Walker spacetime use the 
Arnowitt-Deser-Misner parametrization and consider only the linear scalar perturbation in the gauge in which the 
inflaton perturbations are set to zero Sip — 21]. The metric to linear order is then given by 

goo = -N^ + 9.jN'N^ gQ^ = g,.jN^ g,, ^ a\l + 211)5,, (31 

5™ = -^ r = —, g^^^—^h, (32) 

and 

A^ = l + 5, iV* = -^ + e5,V-"7^. (33) 

ti ati 

Here the derivative with respect to comoving coordinate x^ is indicated as &{. The relevance of the quantity TZ is 
that it represents the combination of adiabatic scalar metric perturbations that is time-independent for Fourier modes 
outside of the horizon, in the limit k/aH — )> 0. Hence, it provides a convenient description of the initial conditions for 
perturbations from inflation. 

In the interaction picture, the free field TZ can be expanded in annihilation and creation operators a(k) and a(k) ' 
that satisfy [a(k),at(k')] = (27r)3,53(k - k'): 

7e(x,77)=/^ 7^fc(77)e"^-a(k)+H.c. (34) 



. = grjN^ 
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0. N^ 
3 = N^ 


g^^ a2 ^'^" 


N' = 


aH 



, (27r)3 . 

and its time evolution is described by the Mukhanov-Sasaki equation |22l.l23j: 

n'l + 2aH{l + 5 + e)7^; + k'^Uu = (35) 

where 5 — H/2HH and e = —H/H^. In the exponential model for the inflaton potential under consideration here, 
5 — —e and Eq. psp can be solved exactly in the case of constant e. The solution with the correct asymptotic behavior 
inside the horizon is; 

7^.(,) = -.{1 - e)e-/^+-/^-|M. (^ff(i)(_fc,) (36) 

= TZk*Uy{-kri) (37) 

with V = 3/2 + e/(l — e), which satisfies the Wronskian condition TZkTZ'^ — TZ'^JZI. ~ i{H/Lp'Y . The constant coefficient 
TZk* is the value of the mode function outside the horizon, TZkir]) = 7?.fe*('7) for —kr] ^ 1, 






(38) 



Since the ?7-dependence above cancels out, H{r]j){—kr]j)^~'' = H{r]^){—kr]t,)^~'' , in which case the mode function 
for Fourier modes outside the horizon may be recast in terms of some other reference time or scale indicated by the 
asterisk: 
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(39) 



where we have also used ?/* = — 1/(1 — e)a*-ff*. (A pedagogical discussion of this model, including the background 
evolution equations and the behavior of the gauge invariant TZ can be found in 2J], sections 4.2 and 10.1, although 
our conventions are slightly different.) The two-point correlation function is used to define the power spectrum P-jz: 



(7^(x,77)7^(y,77)) = 



_fk_ 

(2^ 



e*(^-alp^(fc), p^(fc) ^ |7^,|2 



The power spectrum for modes outside the horizon is 



A^{fc) ^ ^PTiik) = a 






a^^Ht, 



ns-l 



(40) 



(41) 



where C^ = ^^''-^{l - tf^'-^ViyY j-n^ and ns = 4 - 2j/ = 1 - 2e/(l - e). The 7-year WMAP data [251 constrain the 
parameter values ng — 0.967 ±0.014 and C„Hl/eMpy = 2.43(±0.091) x 10~^ at an inverse-length scale a.fH^{^ fco) = 
0.002 Mpc"\ so that e = 0.0162 ± 0.0067 and H^ = 1.12(±0.23) x 10"^ i\/pi. 

III. CORRELATION OF MAGNETIC FIELDS AND CURVATURE PERTURBATIONS 



We now evaluate the cross correlation between the primordial magnetic field and the scalar curvature perturbation, 
TZ, making use of the in-in formalism [2l|, [2a, [23. The expectation value of a cosmological observable O in the state 
of the universe is given in perturbation theory by the standard rules of the in-in formalism: 



{Onit)) 



Texp 



dt' H,r^t{t')] 0'-^'{t)Texp 



dt"H„u{t") 



(42) 



where Oh is the Heisenberg picture operator, Hint is the interaction Hamiltonian and O*"* is the field in the interaction 
picture. (The time integrals are finite in perturbation theory after an appropriate Wick rotation or an analytic 
continuation into the complex plane that introduces an effective small positive imaginary part ie' to time, which we 
indicate with the lower limit 00+.) 

The interaction Hamiltonian is obtained by expanding the action to linear order in perturbations. We work in the 
so-called Maldacena gauge where Sip = 0, which is convenient because perturbations of W vanish. Then, expanding 
the interaction term in Eq. ([1]) to linear order in TZ, 



^int — . I Ct ^ 






tz + ^-]a',^~tzfI 

an , 



^A'^F,,+ed.{W-'TZ')A',F, 



3^ y 



(43) 
(44) 
(45) 



where we have used Eqs. (I5T|) -(|55 )) . Putting these pieces together, and neglecting the term of order e, the interaction 
Hamiltonian in conformal time is given by 



Hi, 



d^x W[f^) 



TZ 



TV_ 
aH 
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If^. 

2 '^ 



2(^^,,p 



(46) 



This Hamiltonian differs from that used in the toy model of our earlier work. In particular, if we eliminate the TZ' 
and diTZ terms and substitute (5(^ for TZ, then Eq. (P5| would resemble the earlier Hamiltonian. 



The cross correlation is 



(7^(x,ry)B,(y,ry)B'(z,J?)) =-2Im / drj,{H,^tirj,)TZi^,rj)Bkiy,v)B''{z,r])) . 

J -00+ 



(47) 



In terms of Fourier-space three point function P3 , 



(7^(x,7y)i3,(y,77)i3'(z,77))=y 



n 






-^ki.x-»k..y-,k3.z (2^)3j3(j^^ ^ 1^^ ^ k3)P3(fcl, fc2, ks; 7^) , (48) 



with 



P3(fci,fc2,fc3;77) ^ -2Im 






drjiX 



(49) 



where the external lines are taken outside the horizon. Then, from the Hamiltonian, the structure of I is found to be 

KM) 



I = w{m) 



T^kiim) 
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a{rji)H{r]i 



{2k2k3fivl^ imyk., im) - ^1^3(1 + tJ'^)vk2 im)vk3 M) 



zKi (?/i) (^(^2, fcs, f^)vk2 iviWk^ im) + ^{h, k2,f^)v'k^ {mW:> im)) 



(50) 



where we have introduced /i — cos6 — k2 ■ ks and '0{x,y,iJ,) = 1 + /x^ + 2/i(a;/j/). Details of this calculation can be 
found in the Appendix. With the aid of Eqs. ([57]) & ([TB)) . one can see that P3 will contain a factor |7?.fci*p|wfc2*p|'yfe3*p 
and a sum of the real parts of time integrals of products of Hankel functions and their derivatives. Specifically, 

P3{ki,k2,k3;T]) = |7^fclH.^Wfe2*^'^fc3*P 



2n/(l-c) 



a4(?7) 
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(51) 
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CtTl 



(53) 
(54) 



where X2 = k2/ki, X3 = k3/ki, r = -fciry and ^36 = K3a{k2 ^ fcs). 

For general slow-roll parameter e, these integrals cannot be evaluated in an analytical closed form, and require 
numerical evaluation. However, the calculation is transparent for the n = 2 case neglecting the slow-roll parameter e 
in the order of the Hankel functions, i.e. we take v = 3/2 and a = 5/2 inside the mode functions. In the appendix 
we list all integrals and their explicit form in this case. 

The amplitude P3 at the end of inflation, in the case ub — and e <C 1, is 



P3{ki,k2,k3,m) = \nk,.\^\vk2*\^\vk,.\''WiTf/a\r^i) 

X {flKi + (1 + fI^)K2 + 19(X2, X3, fi)K3a + ^(^3, X2,lj)K3b) 



Ki 



4 k^k^ 

1-4 \-k2k3{uj'^ + kiuj ~ 2kj) + cjIoj^ - kiuj'^ + 3k?uj ~ fc?) 

9 k^uj-^ 



21,2 r 



K2 = 2- 



- '^■2 "-3 



- (7 + In[-W77/]) + ^^3^ {klkliu'' + fciw - 2kl) 



fci 
-ik2k3Uj{uj^ ~ kiuj^ + ikfio - kf) - 3kiuj^{uj^ ~ Akiu^ - 2kluj + kf))] 



(55) 
(56) 

(57) 



K3a = --r|-4[3cj4-3w3(fci+A:2)+cj2(3fc2 + fc2)+wA:2(3fc2 + fcifc2 + fc2) + 2fcifc|(fci + fc2)] (58) 
and i^36 — K3a(x2 O 0:3), where ui — k\+k2^ k3 and 7 is the Euler-Mascheroni constant. 



We note that, although the individual Fourier modes contributing to the three-point function are constant outside 
the horizon, the correlation continues to grow logarithmically. In our case, the contribution denoted K2 introduces 
such a correction surprisingly at tree level. The term K2 is the analogue to the log{k / aH) contribution found in our 
previous work for the cross correlation of a spectator scalar field with the magnetic field energy density — see I2 for 
n — 2 studied in [2,] . Both of these results are rather puzzling since one would expect log corrections to appear only 
from loops. 

It is appropriate to pause here and note that Eqs. (|51II54I) for the amplitude P3 will be valid even in the absence 
of an amplification mechanism, i.e. even when W{ri) = 1. The result in that case can be obtained by simply setting 
n = and Wj = 1. Carrying out this calculation, we find 

Ki = 2K2 - 4^/(w - fci)(cj + 2ki) X (I + 0{t)) (59) 

Ksa ^ -2^{u + 2ki)x(l + 0{T)), K:ib = K3a{k2^k3). (60) 

However, the sum yields P3 = to the order of approximation of our calculation. Since t <C 1 for modes outside the 
horizon at the end of inflation, any non-zero cross-correlation is exceedingly small in the absence of an amplification 
mechanism. 

IV. BEHAVIOR OF THE CROSS-CORRELATION 

We proceed to analyze the cross-correlation between the primordial magnetic field and the curvature perturbation, 
to determine if there is any imprint of unique signature to indicate the amplification mechanism. 

A. Real-Space Cross-Correlation Coefficient 

The cross-correlation amplitude, evaluated in the coincidence limit, can be determined as follows. Starting from 
Eq. (|55|) . we evaluate the ki integral to eliminate the delta function. The remaining integrand depends only on the 
magnitudes fc2, fca, and 0, the angle between the two vectors: 

{TZB^) = -i^ f kldk2kldk3d{cos9)Ps{kuk2,k3) (61) 

where fci = (fc| + fc| -I- 2fc2fc3 cos 6*)^/^. However, we can replace the 9 integral by fci, whereby 

{TZB^) = -i^ /fcadfcaMfca / k^dk^P3{k^,k2,k3). 



(62) 



Since the integrand is invariant under the exchange of ^2 and fcs, we can replace P3 — >■ 2P-i9{k2 — fca) and remove 
the absolute- value sign from the lower limit of integration. We implement cutoffs at both large and small fc, for the 
ultraviolet and infrared divergences that arise in both the curvature perturbation and magnetic field spectra. The 
cross correlation evaluates to 

where Ni is the number of e-foldings of inflation, r — fcmax/fcmin, and fcmax and fcmin are upper and lower bounds on the 
run of wavevectors. In practice, we expect to link the minimum wavevector with the Hubble scale, fcmin — 27riJo, and 
the maximum wavevector with some galactic scale, fcmax — Stt/A where A ^kpc. Since |fc?7/| ^ 1, we have discarded 
subdominant terms from the above results. The dimensionless cross-correlation coefficient X'jig2 , formed from the 
ratio of the cross correlation with the root-mean-square amplitudes of the curvature perturbation and magnetic fields 
gives 



/T? P'^\ R / FT \ 2 

^ns^ - (^5yT7I(^ = ^ [^J Vl^(^^ - 3 logo 



(64) 



Considering a sufficiently wide range of scales, e.g. r > 10*, and using Nj ~ 70 then X-ji^'^ ~ 900 x (Hj/y/eMpi), 
which is nearly three orders of magnitude larger than a naive expectation for the amplitude. 
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B. Discretized Fourier-Space Cross-Correlation Coefficient 

We now evaluate the triangle-shape dependence of the full three-point correlation function in discrete Fourier space. 
As discussed in Sec. IV C of [7j, this leads to a visualization of the cross correlation that has a clearer imprint of the 
amplification mechanism (cf. Fig. 1 of ibid.). We evaluate a ratio of the form, 

P3{ki,k2,h) ,ggx 



VP-Riki)PBik2)PB{k3y 

to normalize the cross-correlation power spectrum. Since this ratio is not dimensionless, however, we convert the 
continuous Fourier transforms into discretized Fourier transforms. 



d^k 



^E' (66) 



(27r)3 V 
and likewise replace the Dirac delta function with a Kronecker delta, 

(2nfS{ki + k2)^Vdn,,n,- (67) 

We presume a maximum length, L, so that the volume is V^ = L"^ and mode numbers are ki = 2Trni/L. The 
curvature-perturbation and magnetic-field power spectra are now 

(7^2) = ^ e"-(^-^^)/^PK, (68) 

n 

Pn = V'^Pn, (69) 

(b2) ^ ^e-(-^-^^/^<53„s,PB, (70) 

n 

Pb = V-'Pb, (71) 

so that P-ji is dimensionless and Pb has units of (energy)"*. The three-point function becomes 

(TZB^) = y^ gi{ni-x+n2-y+n3-z)/Lp 

ni+n2+n3=0 

h - V-^Pi, (72) 

where P3 has units of (energy)"*. We can now build a dimensionless cross-correlation coefficient, 

C = ^3(ni,n.,n3) _ ^^3^ 



Pn{ni)PB{n2)PBin3) 

where rii for i = 1,2,3 are the magnitudes of vectors Hi that form a closed triangle. 

For isosceles triangles with n2 = n^, the correlation C has two interesting limits. First, for a squeezed triangle, 
with 1 < ni <g; n2 or 6* = TT, then 

Hence, the curvature perturbation and magnetic fields are anti-correlated for squeezed triangles. Second, for a flattened 
triangle with n2 = ni/2 or 9 — 0, then 

C(0 - 0) ^ -^ (^) X 24 (l„(2*u) + . - I) (75) 

where 6 = —2'Krij/L <C 1. Since 6ni ^ 10^^^ for horizon-sized modes, then it is clear that the cross correlation for 
flattened triangle configuration is not only positive but much greater than for the squeezed triangle. A plot showing 
the angular dependence of the cross correlation is shown in Fig. [T] 
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FIG. 1: The ratio C(6')/|C(7r) 
The flattened triangle, at cos 6 



is shown as function of cos 6. We have set 5 
= 1, yields a cross-correlation amplitude C{9 - 



= 2Tvni\rii/L\ ~ 10" 
0) ~ 1500|C(7r)|. 



for horizon-sized modes. 



To show the full Fourier-space triangle dependence of the cross-correlation, we define the quantity 



nsj \C{n)\ 



(76) 



and introduce the variables 0:23 = ^2/71^ and 0:13 = ni/n^, where < X23 < 1 and 1 — 0:23 < 2:13 < 1 + 2^23 covers the 
full set of triangles. As seen in Fig. [2l the maximum value of R occurs for the flattened triangles, corresponding to 
the line X13 = 1 + X23, along which = 0. Squeezed triangles, where 9 = n, are located along 2:13 = 1 — 2:23. 



CONCLUSIONS AND OUTLOOK 



Understanding the origin of galactic magnetism is one key aspect of the evolution and structure of galaxies, and has 
attracted much interest. Single field slow-roll inflation can in general accommodate a natural extension of electrody- 
namics that leaves a random cosmic magnetic field outside the horizon at the end of inflation that is large enough to 
seed the galactic dynamo. In this work we studied the dominant primordial non-gaussian contribution to the statistics 
of this field. This comes from the cross correlation of metric scalar perturbations and magnetic field energy density, 
{TZB'^). We found that in power-law inflation when electrodynamics has a dilaton coupling of the inflaton background, 
there exists a spectrum AfiBB of order 10~^, which is found to peak for flattened triangles in Fourier space. 

We expect two classes of experiments may be sensitive to this primordial non-gaussianity: deep-field all-sky Fara day 
rotation measures and polarization detection of the CMB. For the first class of experiments, examples are SKA |29l - [3l| 
and LOFAR [2^,[34| which have projected sensitives close to 10~^° G. In the SKA experiment, the relevant effect arises 
from the fact that 7?, is directly proportional to the dark matter contrast, hence {TZB'^) constitutes the primordial 
cross correlation of matter density and cosmic magnetic fields. Preliminary searches and studies of this correlation 
include the Refs. ^^. 

In the CMB, cosmic magnetic fields affect most prominently the polarization through the Faraday effect. The 
detectability and bounds on cosmic magnetic field power spectrum in the CMB have been extensively studied in the 
literature: a partial list of recent work relevant to Mpc wavelength fields is Refs. [34-39] . Noting that the temperature 
fiuctuation of the CMB is largely determined by 7?. whereas the B-fields produce a frequency-dependent rotation of 
the polarization along the line of sight, then the cross correlation (TZB^) introduces a primordial non-gaussianity in 
the temperature-polarization-polarization correlation. The feasibility of detecting such a signal remains to be studied. 
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R 



1500 



1000 



500 




FIG. 2: The quantity R, defined in the text as the ratio of the Fourier-space cross-correlation coefficient to that of the squeezed 
triangle, times a factor a::23, is shown as a function of the triangle side lengths. We have set 3 = 2Tmi\rii/L\ ~ 10~ for 
horizon-sized modes. The flattened triangle, at cos9 — 1, yields a cross-correlation amplitude C{9 = 0) ~ 1500 |C(7r)|. 



Appendix A: Calculation of X 



We simplify our calculation by constructing the following differential operators. First, to obtain the magnetic field 
correlation from the vector potential correlation, we introduce 



q4 \^ dy^dz^ dy"^dz'- 



such that B{y) ■ B{z) = O^^Ai{y)Am{z). Second, to obtain the interaction Hamiltonian density, 
(^1) 2^ a{fi,)H{r,i)dm'dmdr,3 









1 



(1 



1 



2' a{r]i)H{r]i) drji 8x2 dx\ 
1 d 



k a^l \ I 



2a{rii)H{rii) dx^ 



a?72 9^3 drji 8x2 



such that 



^™* = X2^^x, (^(ffi) + ^Im + ^Im) 7^(xl,?7l)A,(x2, 772)^5 (x3, 773) 

where Hint — J (Pxia^{'r]i)'Hint- Using these definitions, we can rewrite the cross correlation as 
(7^(x,7,)B,(y,7,)B*(z,77)) = -2Im f dr^i (//,„, (771 )7^(x, 77)5,. (y, 77)B'=(z, 77)) 

J-OC + 



(Al) 

(A2) 
(A3) 

(A4) 

(A5) 
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= -2Im 



00 + 



The expectation value in the above expression breaks into (7?.(xi, ?7i)7?.(x, 77)) and {AiAjAiAm)^ where 



and 



(7^(xl,r,l)7^(x,77)) 



{A{'^2,V2)Aj{-X3,r]3)Ai{y,ri)A„,{z,r])) 



iki-(xi-x) 



+ VkAV3)vkAV2)e''''■^'''-y^+''''■^'''-^^P,l{k2)P^n^Ch) 



(A6) 



(A7) 



(A8) 



with Pij (k) = Sij — {k)i (k)j . In the above expressions and below, a superscript asterisk indicates the complex conjugate 
of a number, not to be confused with a subscript indicating a quantity outside the horizon. 

We can now begin to evaluate the cross correlation. The operator Ob acting on the vector potential four-point 
function gives 

d^k2 d^ks 



o'£'{A,{x2,m)Ai^3,v3)My,v)Amiz,ri)) 



X 'vkAv2)vkMy'''-^'''-''^^''''-^''''^HChUkh - ^^s,J) 

+ «fe.(%)«fe3(^2)e^'^-(''-^-^)+"'-("^-^n(fc2).(^3), - /Aj)] (A9) 

where fJ. — k2 ■ ks — cos 9. Next, including the coupling function W and applying the Oh operators, 

3 

4^1,, E^^H„)'^B"W^('7l)(^(xl,^l)7^(x,,7))(A,(x2,r;2)A,(x3,r;3)A(y,77M™(z,77)) 



»)2,3-s-'Jl n=l 



d^ki d^k2 d^k 



(27r)3 (27r)3 (27r)3 



^ jki (xi -x)+ik2 ■ (xi -y )+ik3 ■ (xi -z) 



T^liiv)v*kJ,v)v*k,{v)^ 



(AlO) 



where 



W{rJl)l2^lk2k 

{l + ^I^)kikl 
1 



Rk^ivi 



RkAm) 



KM) 

aiHi 

KM) 

aiHi 



^LivMkM) 



vk2{'ni)vkM) 



kjklRk^ (?7i) [vl^ {m)vk3 (??i)^(fc3, ^2, A*) + «fc2 (m)"fe3 im)'&if^2, fcs, m)] 



and we define 'i?(a;, y, /i) = 1 + /i^ + 2^{x/y). The cross correlation is then given by 

^g.k,(x,-x)+«k..(x,-y)+.k3.(x,-z)^*^(^)^*^(^)^*^(^)^^ 



(All) 



(A12) 



Rearranging the order of integration, and evaluating the xi integral, we may further simplify this expression to 

(2^ (2^ (2^'"*'"""'"''"*"'^^^)''^'(''^ + ""^ + k3)P3(fci, k2,ks; rj){Al3) 

(TZl (ri)vl (n)vl iri) /"' 1 

P:i{k^MM■,v) - -2Im | ^ „;.^' ' / ^^71! ^ • (AM) 

Hence we arrive at Eq. (|49p . 



«n^) 



00 + 



14 
Appendix B: Integrals 

The integrals required to evaluate the terms Ki in Eqs. (|52ll54p are given as follows: 



Jla = Im <^ / dTiTi '^ -7—U5{x2Ti)-—U5{x3Ti)u3{ti) 



^2_ (,^3 ^ (fc^ ^ fc3)(A;2fc3 - w') + ^^(^2^ + fcl)) X (1 + 0{t)) (B1) 



Jib = Im|y^ dnn '^W|(x2ri)— U5(a;3Ti) (^n— j «3 (ti)| 



^^''s {u;^+uj{k2 + k3) + 2k2h)x{l + OiT)) (B2) 



J2a = Im < dTlTl^^U5{x2Tl)u5{x3Tl)u3{Tl)[ 



= -- (7 + ln[ra;/fci]) + —3-^ (fc^fc2(fci + cj) - 3w/c2A:3(^^ - huj + kl) + 3kiLu^{2ki - w)) x (1 + 0(r))(B3) 

J2b = Im I dTlTl~^U5{x2Tl)u5{x3Tl) fri— )u3(ti) 

= - (7 + ln[Tw/fci]) + -^ {9lj^ - 3uj^{kl + 4fc2fc3 + k^) + Swfcifca/ca - ^k^kl) x (1 + ©(t)) (B4) 

J^a = Im I dTlTl~^U5{x2Tl) (ti—U5{x3Ti)J U3{ti) 

1.2 

3 {3Lu^-3ik2 + k3)uj^ + ikl+3kl)Lu^ + k2{kl + k2k3+3kl)Lu + 2klk3{k2 + k3)) x{1 + 0{t)) (B5) 



9fc3a;3 
J3{, ^ Im < dTiri"''u5(a;3Ti) ( Ti— U5(a;2Ti) j U3(ti) ^ 



= — f^ (3^4 - 3(fc2 + fcs)^^ + (fcg + 3fc^)w2 ^ ^^(^2 ^ ^^^^ ^ 3^2)^ ^ 2fc2fc2(fc2 + ^3)) X (1 + 0{t)) (B6) 

Hence, iCi = -4fcifc2fc3(Jia + Jifc), A'2 == 2{klkl/ki){J2a + J26), and i4r3Q^b = 2(fc|fc|/fci) J3a^b. In all cases we have 
assumed r <C 1. 
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